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SUMMARY 

A droplet camera developed at the NACA Lewis laboratory was used 
to obtain drop-size distribution and drop-velocity data for isooctane 
injected from a sinple orifice directly into a turbulent air stream. 
From these data and wet -bulb tenperature data, vaporization rates and 
drag coefficients were calculated for isooctane drops accelerating and 
evaporating in streams having velocities of 140 and 180 feet per second. 
Spray vaporization rates based on the mean drop diameter D 2 Q were com- 
pared with single-drop vaporization rates by use of the heat-transfer 
equation 
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where DgQ is defined as the diameter of a drop having an area equal to 
the ratio of total area to total number of drops formed in the spray, 6 
is the time, kg is the thermal conductivity of the air. At is the dif- 
ference between the air temperature and drop-surface (wet-biilb) temper- 
ature, is the density of the liquid, is the latent heat of 

vaporization, and Nu 2 q is the heat-transfer Nusselt number based on 
the mean drop diameter D 2 q. For these tests, the change in sensible 
heat of the liquid was negligible con^iared with latent heat requirements . 
The weight percent of the spray evaporated at distances of 1 to 18 in- 
ches downstream of the injector was determined and found to be in agree- 
ment with data obtained in a previous investigation by means of a 
sampling -probe technique. 


Niikiyama-Tanasawa and log-probability drop -size -distribution func- 
tions gave good agreement with experimental drop-size-distribution data 
in the determination of mean drop size. The Sauter mean diameter for the 
spray, predicted by the Nukiyama-Tanasawa empirical expression for a 
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different type of nozzle^ agreed with that found for the 180-feet-per- 
second air stream. The maxim-urn drop diameter for the isooctane sprays 
was calculated from a balance of aerodynamic press"ure and surface ten- 
sion forces under the condition of maxim-um acceleration of the sprays. 
These values compared well with the maximum drop diameter observed from 
photomicrographs of the sprays. 


INTRODUCTION 

When a liquid fuel is injected into the preignition zone of a ram 
jet or a turbojet engine or into an afterburner^ the liquid is broken 
up into a cloud of droplets which are then accelerated to stream veloc- 
ity. While accelerating^ they evaporate at a rate determined by: (l) 

the air-stream velocity^ temperature^ and static pressure, (2) the 
droplet velocity, temperature, and diameter, and (3) the physical prop- 
erties of the liquid and vapor. Thus, the vaporization rate of sprays 
in air streams of known temperature, press-ure, and velocity may be de- 
termined by use of heat-transfer equations when data on drop-size dis- 
tribution, drop acceleration, and drop-surface temperature are available. 

Several investigators (e.g., refs. 1 and 2) have obtained drop- 
size-distribution data for the atomization of liquids in air streams. 
However, the problem of correlating drop-size -distribution parameters 
with physical characteristics of the atomization process is still un- 
solved. Data have been obtained (refs. 3 and 4) for the weight percent 
of fuel spray evaporated in air streams at given distances downstream 
of the injector by means of a sampling technique. Since this method 
did not give data on drop size and drop velocity, a correlation with 
characteristics of momentum, mass, and heat transfer could not be made. 
Correlations of this type have been made for single drops evaporating 
in air streams (refs. 5 and 6). However, the problem of relating 
vaporization rates of single drops to that of sprays has not been solved 
because of the lack of spray vaporization-rate data. 

In reference 7, theoretical equations are given for predicting the 
vaporization rate of fuel spray based on a single mean-diameter droplet. 
For simplicity, an evaporation constant, assumed to be independent of 
drop size and velocity, was used. An expression valid only for either 
exceedingly fine sprays or sprays which have accelerated to stream ve- 
locity was obtained. If the doubtful assumption is made that evapora- 
tion prior to obtaining either of these conditions is negligible, the 
equation of reference 7 may be applied directly to the evaporation of 
fuel sprays in jet engines. 

In this investigation, data on drop-size distribution, drop veloc- 
ity , and drop-s-urf ace (wet-bulb) temperature were obtained for isooctane 
sprays in turbulent air streams having velocities of 140 and 180 feet 
per second. This was accomplished by using a droplet camera developed 
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at the NACA Lewis laboratory which gave photomicrographs of the sprays 
and also drop-velocity data. Drop-size-distribution data were ob- 
tained from photomicrographs of the spray; and the wet -bulb temperature 
for isooctane ^ under the conditions of this investigation, was obtained 
by the experimental technique described in references 5 and 6. From 
these data, it was possible to calculate the following spray parameters: 
(l) mean-drop sizes, (2) vaporization rates, (3) drag coefficients, and 
( 4 ) weight percent evaporated. 


SYMBOLS 


A 

A’ 


a 


b 


The following symbols are used in this report ; 
siarface area of drop, sq cm 
frontal area of drop, sq cm 
acceleration of drop, ft/sec^ 
constant 

molecular mass diffusivity, g/(cm)(sec) 
drag coefficient 


c 

D 

AD 

D 


^20 •> ^30 > 
% 1 ^ %2 


E 

F 

G 

g 


root-mean-square velocity of air molecules, cm/sec 

drop diameter , cm 

drop -diameter size range, cm 

constant 

ZnD’^ 

mean drop diameter defined by the expression D^ ^ 

2nD^ 

for example y D 20 is the diameter of a single drop hav- 
ing an area equal to the total area to nimiber ratio for 
all the drops formed in a spray, cm 

portion of total spray which has evaporated, percent 

total force of air stream, dynes 

momentum 

acceleration due to gravity, 980 cm/sec^ 
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Hv 

h 

k 

L 

I 

M 

m 

dm/ d0 

Am 

Nu 

Nu20 

n 

Pa 

PL 

Q 

q 

R 

Re 

T 

t 

At 

U 


latent heat of vaporization, g-cal/g 
heat-transfer coefficient, g-cal/(sec) (sq cm)(^C) 
thermal conductivity, g-cal/(sec) (sq cm)(^C)/cm) 
distance from mirror to film, ft 
mean free path of air molecules, cm 
magnification factor for lens system, 21 power 
mass of liquid, g 
vaporization rate, g/sec 
mass of liquid vaporized, g 
heat-transfer Nusselt numher, hO/k 

heat-transfer Russelt number based on mean drop diameter 

^20^ ^20 A 

number of drops in given size range 
air-stream static pressure, in. Hg abs 

fuel injection pressure drop, Ib/sq in. 

volumetric flow rate, cc/sec 
constant 

volume fraction of drops having diameter > D 
Reynolds number based on drop diameter, DAup/|i 
temperature , ^R 
temperature, 

difference between air temperature and surface temperature 
of drop, tg - t^, C 

velocity, ft/sec 
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Au 

relative velocity of air with respect to drop, cm/sec 

V 

volume fraction of drops having diameter < D 

We^ 

critical Weber nimnber, Pg^u|Dj^/2a 

X 

distance downstream from the injector orifice, in. 

w y 

C\D 

CD 

w 5 

IniD/B) 

constant 

e 

vaporization time, sec 

X 

evaporation constant, 4kg At/pjE^, sq cm/sec 

M- 

fluid viscosity, poises 

P 

fluid density, g/cu cm 

a 

surface tension, dynes/cm 

0) 

mirror speed, rpm 


Subscripts : 


a 

air 

d 

drop 

f 

final condition 

g 

gas 

i 

initial condition 

L 

liquid 

m 

tnaximuni 

s 

air stream 

V 

vapor 
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APPARATUS AND PROCEDURE 

The vaporization rates of isooctane sprays evaporating and acceler- 
ating in turbulent air streams were determined experimentally with the 
apparatus shown in figure 1. Air at 82°±3° F and 16±1 percent relative 
humidity was supplied from the central laboratory system at 40 pounds 
pej square inch gage. The air was metered with a variable-area orifice 
and exhaus'' d to the altitude exhaust system. In order to minimize 
approach- fa bream turbulence^ a 30 by 28 mesh, 0.013- inch-diameter monel 
wire screen was placed 5 feet 4 inches upstream of the inlet to the 
transparent 8- inch- inside-diameter plastic test section. Two trans- 
parent 2- by 2- by 3/8-inch optically flat plastic windows were installed 
in -.he test section for photographing the spray. Air-stream pressure 
ana velocity in the test section were controlled by upstream and down- 
stream valves . 

Isooctane fuel (2, 2, 4-trimethylpentane, A.S.T.M. specifications) 
was injected by pressurized nitrogen into the test section through a 
1^4- inch-diameter Inconel tube sealed at the end. A 0.041- inch-diameter 
orifice was drilled and reamed to size 1 inch from the sealed end in 
order to obtain a smooth surface. The orifice was pointed directly into 
the air stream and positioned on the test-section center line (position 
4, fig. l) and also moved vertically to positions 1, 2, 6, and 7 in fig- 
ure 1 so that photomicrographs could be obtained over a vertical trans- 
verse of the spray cross section. 

Drop-size and drop-velocity data were obtained by photographing the 

spray at distances of 1, 5^, 14, and 18 inches downstream of the fuel 

injector center line. A camera unit developed at the NACA Lewis lab- 
oratory (ref. 8) was positioned at the 2-inch-square plastic windows, 
as shown in figure 1, to obtain photomicrographs of the isooctane spray 
in natural suspension in the high-velocity air stream. The object plane 
was centered along the test section center line (position 4, fig. l) and 
also moved horizontally to obtain a horizontal traverse of the spray 
cross section (positions 3 and 5, fig. l) . Thus, by adjusting the cam- 
era unit in a horizontal direction and the fuel injector in a vertical 
direction, a seven-position traverse was made of the spray cross section 
at each of the four distances downstream of the injector. 

The lens system used in the camera unit had a magnification factor 
of 21, a 1-centimeter-square field area, a field depth of 1 millimeter, 
and resolved droplets as small as 5 microns in diameter . Stopped images 
of the hi^-velocity droplets were obtained on the film, as shown in 
figure 2, by means of the rotating mirror and synchronized flash system 
shown in figure 1. 
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The original droplet camera described in reference 1 was modified 
considerably in this investigation^ since the additional factors of high 
droplet concentrations and varying droplet velocity were considered. 

The three-phase 400-cycle inverter was rewired and a rheostat was con- 
nected in the field circuit so that the speed of the rotating mirror 
could be varied from 3^000 to 11^000 rpm. The pair of 0.05-microfarad 
condensers were replaced with a pair of 1.0-microfarad condensers which 
gave a higher intensity flash when discharged across the two 0.014- by 
0.1-inch magnesium-ribbon electrodes spaced 3/8 inch apart and sand- 
wiched between split cylinders of bakelite insulation. Also, the lens 
system was redesigned and a modified aerial camera with the shutter and 
lens removed was used primarily as a holder for the 9 .5- inch-width roll 
film. 


By using 1.0-microfarad condensers, which gave a very brilliant 
flash, it was possible to test various fuel dyes, light filters, film 
types, and developing solutions for obtaining high-contrast negatives. 
Best results were obtained with a clear fuel (not dyed), a Wratten red 

No. 25 mounted on the light source, and a shellburst panchromatic- 
type film developed with a D-11 solution. The red light, because of its 
relatively long wave length, gave a minimum of light scattering in the 
spray, and thus sharp droplet images were obtained on the high- contrast 
®Nellburst film. The D— 11 solution was a standard developer generally 
used for hi^-contrast work. 

The following procedure was used in obtaining photomicrographs of 
the spray: The air velocity in the test section was set at 140 or 180 

feet per second. An air temperature of 82°±3° F and static pressure of 
29.3±0.3 inches of merciiry were used in both cases. Isooctane metered 
at a flow rate of 100 pounds per hour with a rotameter was injected from 
the contrastream orifice into the air stream with a pressure drop of 55 
pounds per square inch across the orifice. Photomicrographs of the 
spray were then taken over a range of mirror speeds. From the photo- 
micrographs, the correct mirror speed required for obtaining stopped- 
droplet images on the film and also droplet velocity was determined. 
Photomicrographs of the spray were taken at the proper mirror speed and 
analyzed to obtain drop-size-distribution data. An average of seven pic- 
tures were obtained at each of the seven points in the traverse of the 
spray cross section, and all the droplets in focus and stopped on the 
film were measured to obtain at least 1400 droplets for each analysis. 

The same procedure was repeated at distances of 1, 5^, 14, and- 18 inches 

downstream of the injector to obtain vaporization- rate data. 

Measurements of droplet diameters were made from the photomicro- 
graphs with an eyepiece having small-scale subdivisions of 0.10 milli- 
meter which gave droplet diameters accurate within microns. The 
camera was calibrated with a Bausch and Lomb microscope slide having 
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small-scale subdivisions of 10 microns. Since 1400 to 1800 drops were 
measured for each analysis of the spray cross section, a statistical 

error factor of approximately 3^ percent was calculated for each cummu- 

lative distrilDution analysis. 

The wet-hnlh temperature for isooctane droplets evaporating under 
the conditions of this investigation was obtained experimentally. From 
the data on drop diameter^ velocity, and wet -bulb temperature, it was 
possible to calculate vaporization rates and drag coefficients for the 
isooctane sprays. 


EESULTS AM) ANALYSIS 

The following steps were employed in the analysis of the data: 

(1) A spray- vaporization-rate equation based on the mean drop 
diameter D20 "w-as written for the case of simultaneous acceleration 
and evaporation of sprays. 

(2) Drop-size-distribution data, obtained from photomicrographs, 

were used to calculate the mean drop size D20 each spray analysis. 

(3) Drop velocities were calculated and used to determine vaporiza- 
tion times and droplet drag coefficients. 

(4) A comparison was made between the predicted and observed max- 
imum droplet diameter in the spray. 

(5) Spray vaporization rates based on mean drop diameter Dgo 

were compared with values predicted by single-drop vaporization-rate 
equations (refs. 5 and 6). 

(6) The weight percent of the spray evaporated at distances of 1, 
si, 14, and 18 inches downstream of the injector was determined and 
compared with values predicted by an empirical expression (ref. 4). 


Vaporization-Rate Equations 

The instantaneous vaporization rate dm/d0 of a drop may be de- 
termined from the heat -transfer equation 


= hAAt = jtDkgAtNu 


( 1 ) 
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when the change in sensible heat of the lig.uid is negligible compared 
with latent heat requirements- In recent investigations (refs. 5 and 
6 ) , the following expression was obtained for the heat-transfer Nusselt 
number : 



where the product of the Reynolds and Schmidt nixmbers is given as a 
single dimensionless nmber equal to pAuD/bg^y^ and the momentum- 
transfer ratio g2/c^ is the ratio of gravitational to viscous forces. 

Equation (l) may be rewritten by expressing the instantaneous 
vaporization rate in terms of the change in drop area with time and 
solving for the Nusselt number as follows : 


dD^ 

4kgAt d0 


C3) 


where Nu is the heat-transfer Nusselt number for a single drop. 

The change in area of the spray A may be written as the sum of 
(l) the area of the completely evaporated drops, and (2) the change in 
area of the partially evaporated drops : 

D= VXA6 Nu D=Dp^ 

AA = n ^ ^ nD^ + mXA0 / n Nu (4) 

D=0 D= /^0Nu 

where X = During each experimental time increment A0, the 

decrease in spray area due to complete evaporation of small drops was 
negligible (approximately 1 percent) compared with the decrease in area 
due to partial evaporation of the larger drops in the sprays. Thus, 
the above expression was rewritten as 

D=*Dnj 

^ = «X XD n Nu (5) 

D=0 


For this investigation, equation ("2) may be written 

Nu = 2 + 0.39 Re*^-^ 

(since “ 2.7X10"^ g/(cm)(sec), gl/c^ ~ 1.024X10"^^, and 

(ka/ky)'^*^ = 1.64.) By combining equations (4) and (5) and dividing 

each side of the resulting equation by > y n, the expression 

D=0 
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20 


A9 



( 6 ) 


is obtained, where 


D^q * 2nD^^2n» 


It was found experimentally that : 


D=D, 



0.90 i 0.02 


( 7 ) 


Thus, combining this expression with equation (6) and passing to the 
limit yields 


^20 

de 


2 + 0.35 


Peo^Pa 

^^a 


0.6 


PL^v 


Nu20 


( 8 ) 


where the term in brackets represents the Nusselt number based on the 
mean drop diameter PgO’ 


The mean Nusselt number Nu20 niay be determined from equation (S) 
and compared with the single-drop Nusselt number given by equation (2) 
when the spray parameters dD^Q, Au, and -At are known. Since the 

surface temperature of each drop could not be determined experimentally, 
the wet-bulb temperature of isooctane was obtained experimentally and 
assumed equal to the surface temperature of a single drop having a mean 
drop diameter equal to D20* The remaining parameters dDgo and Au 
were determined experimentally. 


Drop-Size-Distribution Analyses 


Determination of mean drop sizes . - Analyses of the drop-size dis- 
tribution for isooctane sprays at distances of 1, S^-, 14, and 18 inches 

downstream of the fuel injection station were made for air-stream veloc- 
ities of 140 and 180 feet per second. Table I shows the drop-size count 
obtained at each of the seven positions in the spray cross-section tra- 
verse for a representative test. A complete traverse of the positions 
1 to 7 and total drop-size count are included in this table. Eight 
total-drop- size counts are recorded in table II, and a typical plot of 
drop-size distribution is shown in figure 3. 
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Three of the most common mathematical expressions for drop-size 
distribution are the Rosin-Rammler ^ Nukiyama-Tanasawa^ and log- 
probability equations. A description of the evaluation and use of these 
expressions is given in reference 9. They may be written as follows: 

Ros in-Rammler ^ 


CKl 

iro 

CD 

Rukiyama-Tanasawa ^ 


log-probability, 

^ e-SV (11) 

iy V” 

Drop-size-distribution data were plotted for each of the eight 
analyses using each of the three distribution equations. Samples of 
these plots are shown in figures 4 to 6. The Nukiyama-Tanasawa and 
log-probability expressions gave good agreement with values of D20 ob- 
tained by direct integration of the data. Results for the test condi- 
tion Us = 140 feet per second and x == 5.5 inches are shown in the 
following table : 


Mathematical expression 

Mean drop diameter, 
microns 

Rosin- Rammle r 

21.1 

N ukiyama - Tana s awa 

41.3 

Log-probability 

40.7 

Direct integration 

41.3 




dD 


D 


(9) 


dD 


120 


^5 -bD 
D e 


Cio) 


The Ros in-Rammler expression gave poor agreement since q in equation 
(9) was found to be approximately equal to 3. When q is equal to 3 
in this expression, number mean diameters such as D20 Q*i’e predicted 
to be equal to zero. This, of course, limits the usefulness of equation 
(9) . A similar limitation of the Ros in-Rammler equation was found in 
reference 9 . 

Comparison with predicted mean drop sizes . - Reference 10 gives 
the following empirical expression for predicting the Sauter mean diam- 
eter D32 produced by a gas atomizing nozzle: 
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D 


32 


Us 'VPl [ ac^j \ «s / 


1.5 


( 12 ) 


For this investigation, Qj^/Qg is small, and the expression reduces to 


- 1920 
®32 - 



10,516 

Us 


(13) 


for isooctane. By using the Nukiyama-Tanasawa size-distribution expres- 
sion, D 32 1 ® related to DgQ as follows (ref. 11 ): 


D 20 = 0.693 D 


32 


which may be substituted into equation (13) to give 

7288 


1>20 = 


U, 


(14) 


Values for the mean drop size D 20 obtained 1 inch downstream of 

the injector for the two air-stream velocity conditions are listed in 
the following table, with values obtained from equation (14) for 
comparison: 


Air-stream 

velocity. 

Us. 

ft/ sec 

Mean drop size, microns 

Experimental 

From eq. (l4) 

180 

140 

40.7 

45.0 

40.5 

52.1 


The agreement is good at the higher air-stream velocity. However, the 
effect of Ug on D 20 appears to be greater for the calculated values 

shown in the preceding table. This may be because of the different 
type of nozzle used in the Nukiyama-Tanasawa atomization studies. It 
is, therefore, doubtful that the Nukiyama-Tanasawa equations can be 
used to predict initial drop size over wide ranges of operating condi- 
tions with the simple orifice injector used in this investigation. 


Analysis of Drop-Velocity Data 


Calculation of drop velocity . - The velocity of droplets which ap- 
peared on the photomicrographs as stopped images was calculated from the 
expression 


Ud = 


4itLcJ0 


(15) 
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Sharply defined droplet images were obtained on the film when mirror - 
speed settings gave values of drop velocities calculated by the preced- 
ing expression within 15 feet per second of the actual drop velocity. 
Since the duration of the magnesium flash was 4 microseconds, a droplet 
travel of 6 microns was permissible and droplet images on the film 
could be measured within i-3 microns. Droplet images were badly smeared 
and lacked sufficient contrast for measuring when mirror-speed settings 
gave drop velocities which were not within ±5 feet per second of the 
drop velocity. 

The proper mirror speed was determined by taking photographs at a 
low mirror speed and gradually increasing the speed until stopped drop- 
let images were obtained. The photographs showed that droplets of all 
sizes (table l) were moving at approximately the same velocity, because 
droplets of all sizes were stopped at a particular value of mirror speed. 
It was also established experimentally that a range of mirror speeds 
could be used in which photomicrographs were obtained showing stopped 
images of drops of all sizes at each of the mirror speeds. The range 
of mirror speeds resulting in definition of the drops was greatest at 
the station 1 inch downstream of the fuel injector; at this location the 
range of mirror speeds giving drop definition was about 20 percent of 
the mean value. Approximately the same drop-size distribution and mean 
drop size were obtained from the highest and the lowest mirror speeds 
giving drop definition. This range of mirror speeds was probably due 
to the liquid jets penetrating several inches into the air stream and 
being broken up into drops in the penetration zone shown in figure 7. 
Thus, droplets formed at the maximum point of penetration had time to 
accelerate to a relatively high velocity, compared with drops formed at 
the orifice, before reaching the camera station. Images of droplets of 
all sizes could be stopped at a particular value of mirror speed with 
the injector at each of the different locations. 

Minimum mirror speeds which gave stopped droplet images were used 
in obtaining photomicrographs for each spray analysis so that the in- 
jector orifice could be considered to be the point of initial formation 
and acceleration of the spray. Drop velocities are plotted against the 
distance downstream from the fuel-injector orifice in figure 8 for the 
droplets starting at zero initial velocity at the injector orifice. For 
the case of solids injected into the 180-foot -per-second air stream, a 
5-micron-diameter sphere would accelerate to approximately stream veloc- 
ity while traveling a distance of 1 inch, whereas 5-micron drops were 
found to have a minimum velocity of 85 feet per second, as shown in fig- 
vre 8. Thus, the acceleration of small drops was found to be consider- 
ably lower than that predicted for small solid spheres. 

Calculation of vaporization time . - Vaporization time was calcu- 
lated by first determining an empirical relation between drop accelera- 
tion and velocity difference, and then integrating the expression in two 
successive steps to obtain an equation giving distance of droplet travel 
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as a function of stream velocity and time. Acceleration was determined 
from a plot of drop velocity squared against distance (fig. 9) and the 
acceleration equation 


u: 


a = 


d,f -%i 
2x 


Thus, the acceleration at any distance x was equal to one-half the 
slope of the curve in figure 9. 

In order to determine a relation between drop acceleration and 
velocity difference, the plot shown in figure 10 was prepared. From 
this plot, the following empirical expression was obtained: 

a = 333 (Us - U(i) (16) 


Since acceleration may be expressed as the change in drop velocity with 
time dU^/de, the preceding expression may be written as 



Integrating between the limits shown gives 

U , = U (1 - 
d s '■ 


(17) 


By substituting dx/d0 for droplet velocity in the preceding ex- 
pression. 




e-3330 


d0 


is obtained. 


Integrating gives the final expression 


X = U3 


1 


g-3330 

333 


(18) 


where the vaporization time Q is an implicit function of the stream 
velocity Ug and the distance x. 
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Calculation of drop drag coefficients . - The following force "bal- 
ance for an evaporating droplet accelerating in an air stream was de- 
rived (appendix A) : 


F = ma + Au dm/d0 (l9) 

The instantaneous aerodynamic pressure force of the air stream accelerat- 
ing the vaporizing drop may he expressed as 

F = I PgA'(Au)2Cp (20) 

The following equation for the drop drag coefficient was obtained 
hy combining equations (19) and (20) and substituting into the expres- 
sion tenns for drop volinrie^ drop frontal area, and vaporization rate 
(given by eq. (12) ) : 


4 Da 


8k AtNu 


Cj) = ~ — 


3p /A^2 Hp DAu 
*^g (Au) v"^g 


( 21 ) 


(For this investigation, p^ = 0.69 g/cu cm, pg = 12X10“'^ g/cu cm, 

kg = 5.9x10-^ g-cal/(cm)(sec)(Oc), At = 25^ C, and = 73 g-cal/g.) 

Drop drag coefficients were calculated from the preceding expression 
and plotted against the drop Reynolds number as shown in figure 11 for 
five velocity differences. This figure shows that the droplet drag 
coefficient for evaporating liquids cannot be correlated directly with 
the Reynolds number to obtain a single curve for all velocity differ- 
ences. Drag-coefficient values for solid spheres in air streams (ref. 
12) are also plotted in figure 11. 


An empirical expression was obtained for the drop drag coefficient 
by combining equations (l6) and (2l) to give 


D 

Cj3 = 444 — ^ + 


p Au 


8k AtNu 
g 

H p EAu 
v*^g 


( 22 ) 


where the constant 444 contains the units reciprocal seconds required 
to balance the equation dimensionally. Further tests on the accelera- 
tion of drops are needed using other liquids and varying air-stream 
pressure and temperature conditions in order to establish a more general 
expression. 


Comparison of predicted and observed maximum droplet diameters . - 
When the aerodynamic forces over an accelerating drop exceed the surface 
tension forces, the drop will shatter. Thus, it is possible to calcu- 
late the maximum drop diameter that will be found in a spray under the 
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condition of maximum acceleration at the injector orifice. Some investi- 
gators (refs. 13 and 14) have studied the break-up of water drops, and 
equated aerodynamic pressure forces and surface tension forces as 
follows : 


Cj) i pg(Au)^A' = rtDa 


where A' and D are the maximum drop frontal area and diameter^ re- 
spectively. Maximum drop diameters, calculated from this expression, 
are compared with observed values in the following table, which also 

shows values of the critical Weber number (We = paU^D /2a) : 

c s m 


Air-stream 

Maximum drop diameter. 

Critical Weber 

velocity. 



number, We^, 

ft/ sec 

microns 

based on ob- 
served 

Calculated 

Observed 

140 

127 

126 

6.7 

180 

108 

112 

9.6 


Thus, agreement was obtained between values of Dm calculated for at- 
omization at the orifice and observed values obtained 1 inch downstream 
of the orifice. 


Fue 1 -Spray -E vaporat ion C alcu lat i ons 

Comparison of fuel-spray and single-droplet vaporization rates . - 
From the drop-size distribution, drop velocity, and wet-bulb data, the 
instantaneous vaporization rate of the spray expressed as the change in 
area of the mean drop size D 20 with time may be obtained by'plotting 

D^ against 6, as shown in figure 12. In order to compare single- 

droplet vaporization-rate equations with experimental spray vaporization 
rates, equations (2) and (3) were used to calculate for a 

single droplet, and the results were plotted in figure 12. Experimental 
and predicted values of D 2 Q agree within approximately i5 percent. 

Thus, the assumption that the mean surface temperature is closely ap- 
proximated by the wet-bulb temperature appears valid in this investiga- 
tion. This validity was possibly a result of low over-all fuel-vapor 
concentration. Although the concentration of drops was high near the 
injector orifice, the velocity difference Au was also very high 
(droplet film thicknesses were small) so the spray evaporated essen- 
tially as a cloud of isolated droplets. Simplified vaporization-rate 
equations are derived in appendix B. 
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Comparison of predicted and experimental weight percent of spray- 
evaporated . - Since the change in spray mass due to the complete evap- 
oration of small drops was negligible compared with the change in mass 
due to the par'tial evaporation of larger drops in the spray, the weight 
percent of the spray evaporated E was calculated from the expression 


E = 100 


(^ 50 )j ~ (Pjo^f 

^^30 


•z 

Values of D30; olotained ly direct integration of the drop-size- 

distrihution data^ are given in table III with the weight percent of 
the spray evaporated. Initial values of D20 were obtained from the 

ciiTve in figure 12 at 6 = 0. Using these values, initial values of 
D30 were determined from the expression 

D30 = 1.10 D20 (23) 

Experimental data showed the ratio D30/D20 was equal to 1.10±0.02, and 
the Nukiyama-Tanasawa expression for mean drop sizes gives 

D30 = 1.13 D20 


Table III also gives values for the weight percent of spray evap- 
orated as calculated from the following empirical expression derived in 
reference 4: 


E 

100 - E 


= 9.35 


lOOC/ 


4.4 


100 / 


vO.80 


^-1.2 0.42 x<^.84 


(24) 


The agreement is good except for the station 1 inch downstream of the 
fuel injector, which represents an extrapolation of equation (24), since 
the sampling probe technique could not be used for distances of less 
than 5 Inches from the fuel injector. 


SUMMAEY OF RESULTS 


A droplet camera was used to obtain drop-size distribution and 
drop-velocity data for isooctane sprays. An analysis of heat-transfer 
and drag data for the sprays evaporating and accelerating in turbulent 
air streams gave the following results: 

1. Mean drop sizes were determined by direct integration of the 
data and compared with values obtained from an analysis of the data using 
the Rosin-Rammler, Nukiyama-Tanasawa, and log-probability expressions 
for drop-size distribution. Agreement was obtained with the Nukiyama- 
Tanasawa and log-probability expressions. 
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2. The mean drop diameter D^q predicted hy the empirical Nukiyama- 

Tanasawa expression for a different type of nozzle agreed with the ex- 
perimental value of D^o obtained at a distance of 1 inch downstream 

of the injector in a 180-foot -per -second velocity air stream. The effect 
of stream velocity on drop size predicted in the Nukiyama-Tanasawa ex- 
pression is different^ however^ from that obtained experimentally. 

3. An empirical expression was derived for the drag coefficients of 
isooctane drops accelerating and evaporating in turbulent air streams. 

4. An expression^ obtained by equating aerodynamic pressure forces 
and surface-tension forces for drops^ was used to calculate the maximum 
drop size existing in the spray for the condition of maximum spray ac- 
celeration. Photomicrographs gave maximum drop diameters which were in 
agreement with this expression. 

5. The vaporization rate of an isooctane spray based on the mean 
drop diameter D 2 q (defined as a drop having an area equal to the total 

area to total number ratio for all the drops in the spray) was found to 
correlate single-drop vaporization rates. The weight percent of the 
spray evaporated downstream of the injector was determined and found 
to agree with values obtained using a different technique. 


CONCLQDING REMARKS 

Spray- vaporization-rate data obtained in this investigation corre- 
late single -droplet vaporization rates. The Reynolds number effect 
during spray acceleration was found to be very important in determining 
the time required to evaporate the spray. Approximately 50 percent by 
weight was evaporated during the acceleration period. Thus, when more 
data on spray acceleration and atomization are available, the expres- 
sions used in this investigation may be tested and extended in the field 
of spray atomization, acceleration, and evaporation. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, July 30, 1954 
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APPENDIX A 


OJ 

CO 

CO 


CHANGE OF MOMENTUM OP EVAPORATING DROP 

A mass m + An is considered to be moving at a velocity as 

shown for condition I: 


Condition I 


Condition II 
(S) -►Ug 



The velocity of the mass m is increased to and the velocity 

of the element An is increased to stream velocity Ug, as shown above 
for condition II. Thus^ the change in momentum is given as 

follows : 


A} = m (Ud + AUd) + Am Ug - (m + An) Ud 


which reduces to 


A} = mAUd + An(Us - Ud) 


Dividing by A0 and passing to the limit yields the instantaneous force 
F as follows: 

T, dG *^^d „ •> dm 

^ = le “ d^ (Us - Ud) ^ 

where dm/dS is the instantaneous vaporization rate. 
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APPENDIX B 


FUEL-SPRAY -EVAPORATION CALCULATIONS 

The time required to evaporate a spray may be divided into two 
periods: (l) the period required to accelerate the spray to approxi- 

mately stream velocity, and (2) the remaining period of evaporation in 
which the spray travels at approximately stream velocity. In this in- 
vestigation, the spray acceleration period was studied, and the spray 
vaporization rate (based on the mean drop diameter DgQ) was found to 
agree with equation (s) : 

/D^ 

2 + 0.35\^ 


®20 "^g^^ 

d0 » 


The final period of spray evaporation has been treated in reference 
7 for the spr^y traveling at approximately stream velocity. Thus, a 
complete analysis of fuel spray evaporation may be made by treating 
each period separately and combining the results. 


Evaporation During Spray Acceleration Period 


Determination of change in D^q during spray acceleration period. 


For this investigation. 


kg, g-cal/(sec) (sq cm) (°C/cm) 

At, °C 

PL^ s/cu cm 

Hy, g-cal/g 

and the vaporization-rate equation may be written as 


dD 


20 


d0 


r= 11.65X10-5 


2+0.35 


D^^Aup 
20 a 


0.6 


5.9X10-5 

25 

0.69 

73 


Equation (17) gives Au in terms of Ug and 0. By substituting this 
expression into the vaporization-rate equation. 


dD 


d0 


= 11.65X10-5 


2 + 0.35 


^2oJ^a 


0.6 


-200 e 


(25) 
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is obtained. This expression may be graphically integrated by taking 
small increments of Q and calculating the change in d|o’ 


It was found experimentally that the sprays had accelerated to 
approximately stream velocity in 0.014 second. Althoiigh decreased 

0 

approximately 33 percent in this time^ ^20 decreased only 5 percent. 
Thus^ if 1^20^ assumed approximately constant^ the vaporization rate 
equation may be integrated to give 

Ad|q = 11. 65X10"^ Jo. 028 + 1.75X10"^ 



since, 



, -200 0 
1 - e 

200 


+ C 


where C = 0.06 for the condition of complete spray acceleration. By 
using the following values: 


(d|o)., sq cm 
ft/sec • 
poises . 
gm/cu cm 


2100x10"® 

140 

1.8X10"^ 

12.04x10"^ 


the expression 

ADgo = (326 + 375)X10"® = 701X10'® sq cm 

2 

may be written. Graphical integration of equation (25) gave DgQ = 
700X10"® square centimeter (fig. 12) for Ug = 140 feet per second, and 

0 = 0.014 second. Thus, a simple expression was obtained, since the 
change in the Reynolds number was controlled primarily by the change 
in velocity difference Au and was affected very little by the change 
in D20» 


Calculation of weight percent evaporation during spray acceleration 
period . - The weight percent evaporated E during the spray accelera- 
tion period may be approximated by the expression 


E = 100 X 


(D® ) /(D® ) 

^ 20 ^^^ ^ 20*^1 


1.5 


(26) 


since it was found experimentally that ^zo/^ZO ~ 1 *1010 *02. Values 
of E calculated from this expression and experimental values of E 
are compared in the following table: 
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Air -stream 

Time, 

Distance from 

Evaporated portion of 

velocity. 

0 , 

injector, x. 

spray, E, 

percent 

Us. 

ft/sec 

sec 

in. 

Experimen- 
tal, from 
table III 

Calculated 
from eq. (26) 

140 

0.011 

14 

40.0 

44 

180 

0.011 

18 

51.2 

51 


Evaporation During Period Following Complete Spray Acceleration 

The method of reference 7 may be used to calculate the weight 
percent evaporated after the spray has accelerated to stream velocity. 

In this method^ the Rosin-Rammler expression (eq. (9)) is used. When the 
values of D and q are known, the weight percent unevaporated 100 - E 
may be determined. 

The relation between D and D^q may be derived in the following 

manner: The Nukiyama-Tanasawa equation for mean drop sizes (ref. ll) 

predicts I>3i/l^20 = 1*29 which agreed within 2 percent of experimental 
values of I^3i/^^20 obtained in this investigation. Also, it was shown 

in reference 7 that the Rosin-Rammler expression for mean drop sizes 
gave D = /\/2 D 31 as a close approximation for D when 2<n<4.. Thus, 


may be written. Figure 12 shows D 20 == 43.6 microns at Ug = 140 feet 
per second and x = 5.5 inches (9 = 0.0058 sec_^. By substituting this 
value of DgQ into the preceding expression, D = 79.5 microns is ob- 
tained. This value agrees with figure 4 which shows D = 80 microns. 
Experimental values of D/D 20 agreed with the constant 1.82 within 2 
percent . 

In this study, the constant n in equation (9) was approximately 
3 throughout the spray acceleration period. By using n = 3, and 
calculating D from the expression 


D = 1.82 - 11.65X10-5 1^0.028 + 1.75X10“^ 


^®2oh^sP£ 


0.6\0.5 


The weight percent evaporated E following the acceleration period may 
be determined from the plot of reference 7 of 100 - E against 
(where X = 23.3 10"^ cm^/sec and 9 is the time following the spray 
acceleration period). Thus, the total weight percent evaporated may be 
calculated as the percent evaporated during the spray acceleration period 


3293 


3293 


NACA TN 3265 


23 


plus the weight percent of the remaining fraction evaporated in the 
final evaporation period. 


Application of Method of Reference 7 to Spray Acceleration Period 

If it were assumed that the spray traveled at stream velocity 
(Au = 0 and Nu = 2) in the 140-foot-per-second air stream for the 
distance x = 18 inches (0 = 0.014 sec), the weight percent evaporated 
wo\ild be calculated in the following manner: From the plot of refer- 

ence 7 of weight percent unevaporated spray (lOO - E) against 
the expression 


may be obtained, from which 100 - E = 88 percent, or E = 12 percent 
evaporated. This value of E is approximately 30 percent less than 
the experimental value given in table III. Thus, a very serious error 
would be made if the evaporation of the spray during the acceleration 
period was ass\mied the same as evaporation in a quiescent atmosphere. 
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TABLE I. - DROP-SIZE-DISTRIBUTION DATA FOR SEVEN-POINT TRAVERSE 

OF SPRAY CROSS SECTION 


[Air-stream velocity, Ug, 140 ft/sec; distance from injector 

orifice, x, 18 irQ 


Drop-diameter, D, 
microns 

Number of drops, n 

Traverse position^ 

Total 

for 

first 

trav- 

verse 

Second 

trav- 

verse 

Total 
for two 
trav- 
verses 

Range^ 

Average 

1 

2 

3 

4 

5 

6 

7 

0- 4.8 

2.5 

0 

0 

0 

1 

2 

1 

2 

6 

4 

10 

4.8- 9.5 

7.0 

0 

2 

1 

5 

2 

5 

3 

18 

15 

33 

9.5- 14.3 

12.0 

3 

6 

0 

8 

7 

10 

8 

42 

54 

96 

14.3- 19.0 

17.0 

9 

11 

6 

13 

8 

16 

10 

73 

56 

129 

19.0- 23.8 

21.5 

13 

27 

29 

31 

12 

24 

16 

152 

103 

255 

23.8- 28.5 

26.0 

16 

18 

14 

19 

15 

30 

20 

132 

82 

214 

28.5- 33.3 

31.0 

18 

6 

24 

9 

17 

17 

17 

108 

78 

186 

33.3- 38.0 

36.0 

15 

11 

14 

11 

15 

8 

17 

91 

36 

127 

38.0- 42.8 

40.5 

14 

10 

17 

4 

9 

8 

10 

72 

29 

101 

42.8- 47.5 

45.0 

7 

11 

9 

5 

4 

2 

8 

46 

24 

70 

47.5- 52.3 

50.0 

3 

4 

4 

4 

5 

3 

7 

30 

21 

51 

52.3- 57.0 

55.0 

5 

6 

4 

3 

9 

1 

3 

31 

13 

44 

57.0- 61.8 

59.5 

2 

2 

1 

2 

4 

1 

4 

16 

20 

36 

61.8- 66.5 

64.0 

9 

2 

1 

3 

4 

1 

3 

23 

6 

29 

66.5- 71.3 

69.0 

6 

4 

2 

3 

5 

0 

1 

21 

6 

27 

71.3- 76.0 

74.0 

5 

2 

0 

2 

0 

0 

0 

9 

10 

19 

76.0- 80.8 

78.5 

1 

1 

2 

3 

3 

2 

1 

13 

0 

13 

80.8- 85.5 

83.0 

1 

0 

1 

3 

2 

0 

0 

7 

2 

9 

85.5- 90.3 

88.0 

2 

2 

1 

0 

3 

0 

0 

8 

1 

9 

90.3- 95.0 

93.0 

0 

2 

0 

0 

3 

0 

0 

5 

3 

8 

95.0- 99.8 

97.5 

0 

0 

0 

1 

0 

1 

0 

2 

1 

3 

99.8-104.5 

102.0 

0 

1 

0 

0 

1 

0 

0 

2 

0 

2 

104.5-109.3 

107.0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

1 

109.3-114.0 

112.0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

1 

Total 

129 

130 

130 

130 

130 

130 

130 

909 

564 

1473 


®'AD = 4,75 microns. 
^Fig. 1. 
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TABLE II. - SUMMARY OF DROP-SIZE-DISTRIBUTIOW DATA 


Average 


Air- 

•stream velocity, Us, ft/sec 


drop 

diameter , 

140 

180 

D, 


Distance from injector, x, in. 


microns 









1 

5.5 

14 

18 

1 

5.5 

14 

18 





NumUer of 

drops . 

^ n 



2.5 

8 

12 

10 

10 

16 

20 

23 

26 

7 

24 

37 

36 

33 

30 

49 

65 

54 

12 

51 

76 

78 

96 

66 

66 

102 

111 

17 

76 

120 

148 

129 

86 

103 

183 

188 

21.5 

115 

176 

194 

255 

132 

168 

237 

234 

26 

155 

200 

229 

214 

165 

189 

304 

241 

31 

187 

260 

201 

186 

241 

239 

298 

233 

36 

234 

267 

153 

127 

241 

190 

221 

158 

40.5 

144 

172 

115 

101 

148 

136 

157 

122 

45 

127 

123 

85 

70 

125 

100 

114 

81 

50 

100 

72 

70 

51 

84 

73 

60 

66 

55 

93 

68 

57 

44 

73 

60 

51 

53 

59.5 

74 

60 

37 

36 

57 

38 

36 

28 

64 

53 

40 

30 

29 

44 

25 

25 

25 

69 

44 

36 

20 

27 

25 

20 

17 

12 

74 

33 

23 

17 

19 

18 

15 

11 

12 

78.5 

31 

19 

13 

13 

11 

10 

8 

5 

83 

24 

17 

13 

9 

12 

7 

4 

5 

88 

14 

13 

11 

9 

8 

3 

4 

3 

93 

9 

8 

4 

8 

6 

3 

3 

1 

97.5 

4 

5 

5 

3 

4 

1 

1 

- 

102 

5 

4 

2 

2 

3 

2 

2 

- 

107 

3 

3 

2 

1 

1 

1 

- 

- 

112 

1 

3 

1 

1 

1 

1 

- 

- 

116.5 

1 

2 

- 

- 

- 

- 

- 

- 

121 

0 

2 

- 

- 

- 

- 

- 

- 

126 

1 

1 

- 

- 

- 

- 

- 

- 

Total 

1611 

1819 

1531 

1473 

1597 

1519 

1926 

1658 
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TABLE III. - WEIGHT PERCENT OF SPRAY EVAPORATED 


Distance 

from 

injector, 

in. 

dIqXioIo, 

cc 

Weight percent evaporated 

Experimental 

Calculated from 
equation (23) 

Air-stream velocity, U^, 140 ft/sec 

0 

^1279 

_ _ _ „ 

0 

1 

1182 

7.6 

7.0 

5.5 

952 

25.6 

24.1 

14 

768 

40.0 

41.0 

18 

738 

42.3 

46.5 

Air-stream velocity, Ug, 180 ft/sec 

0 

^1036 


0 

1 

863 

16.7 

8.9 

5.5 

692 

33.2 

28.0 

14 

537 

48.3 

46.0 

18 

506 

51.2 

51.5 


^Calculated from D^o/DgQ = 1.10 where DgQ was 


obtained from fig. 12 at x = 0. 
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Figure 1. - Schematic diagram of test section and camera unit. 
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Figure 2. - Photomicrograph of isooctane spray in 180-foot-per-second air stream 
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Drop diameter, D, microns 

Figure 3. - Droplet-size distritutlon. Air-stream velocity, 140 feet per second; dls 
tance from injector, 5.5 inches- 
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Drop diameter, D, microns 

Figure 4. - Rosin euid Rammler analysis- Air-stream 
velocity, 140 feet per second; distance from injector, 
5.5 inches; D, 80 microns; mean drop diameter 
21.1 microns; q, 3-1. 
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Figure 5. - Nukiyama and Tanasawa analysis. Air-stream velocity, 140 feet per second; 
distance from injector, 5.5 Inches; mean drop diameter D 30 41.3 microns; b, 0.0837. 
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Figure 7 - 


Break-up of liquid fuel jet in turbulent air stream. 
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Figure 9. - Acceleration of Isooctane droplets- 
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Figure 10. - Relation between drop acceleration and 
velocity difference. 


Drag coefficient 



Figure 11« - Relation between drag coefficient and Reynolds number for isooctane 
droplets and solid spheres. 
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Figure 12. - Correlation of spray vaporization rate 
(based on mean drop diameter D 2 q) with single- 
drop vaporization-rate equation. 
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